Nucleation pathways on complex networks 
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Abstract 

Identifying nucleation pathway is important for understanding the kinetics of first-order phase 
transitions in natural systems. In the present work, we study nucleation pathway of the Ising 
model in homogeneous and heterogeneous networks using the forward flux sampling method, and 
find that the nucleation processes represent distinct features along pathways for different network 
topologies. For homogeneous networks, there always exists a dominant nucleating cluster to which 
relatively small clusters are attached gradually to form the critical nucleus. For heterogeneous 
ones, many small isolated nucleating clusters emerge at the early stage of the nucleation process, 
until suddenly they form the critical nucleus through a sharp merging process. By analyzing the 
properties of the nucleating clusters along the pathway, we show that the main reason behind the 
different routes is the heterogeneous character of the underlying networks. 
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I. INTRODUCTION 



Nucleation is a fluctuation-driven process that initiates the decay of a metastable state 



into a more stable one [l 
growth of a new phase [2 



. It is usually involved in first-order phase transitions and along with 
^M|- Many important phenomena in nature, including crystallization 
6|, fractures |3, S], glass formation {q], and protein folding 10 1, to list just a few, are 
associated with nucleation. Despite much attention, many aspects of nucleation processes 
in complex systems are still unclear and deserve more investigation. 

The Ising model is a paradigm for many phenomena in statistical physics. It has also 
been widely used to study the nucleation process. For instance, in two-dimensional lattices, 
Allen et al discovered that shear can enhance the nucleation rate and the rate peaks at an 
intermediate shear rate Sear found that a single impurity may considerably enhance 



the nucleation rate 



12|. Page and Sear reported that the existence of a pore may lead 



to two-stage nucleatiori, and the overall nucleation rate can reach a maximum level at an 
intermediate pore size 13(|. The nucleation pathway of the Ising model in three-dimensional 



lattices has also been studied by Sear and Pan [Ij, |15|. In addition, the Ising model has 
been frequently used to test the validity of classical nucleation theory (CNT) |16|-|22|. Nev- 
ertheless, all these studies are limited to regular lattices in Euclidean space. 



Since many real systems can be properly modeled by network-organized structure 
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25| . it is thus an interesting topic to explore nucleation process in complex networks. Very 



recently, our group have studied nucleation dynamics on scale-free (SF) networks [2&\ and 
modular networks j^^]. In these two papers, we mainly focused on the nucleation rate and 
system size effects. We found that, for SF networks, the nucleation rate decays exponentially 
with network size, while the critical nucleus size increases linearly. For modular networks, as 
the network modularity worsens the nucleation undergoes a transition from a two-step to one- 
step process and the nucleation rate shows a nonmonotonic dependence on the modularity. 
As we know, network topology could play important role in the system's dynamics, involving 
not only the stationary state, but also the detailed pathway. For instance, it was shown that 
network heterogeneity could drastically influence the path to oscillator synchronization 
Nevertheless, how network topology would influence the nucleation pathway has not been 
reported yet. Motivated by this, we will study the different roles of network architectures in 
the formation of nucleating clusters, which can reveal the nucleation pathways of the Ising 



28|. 
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model in the underlying networks. 

Since nucleation is an activated process, it can be extremely slow. Therefore, direct 
simulations can take excessive amounts of time. To overcome this difficulty, in the present 
work, we adopt the forward flux sampling (FFS) [29] approach proposed recently, which is 
efficient and easy to implement to study rare events. We employ Erdos-Renyi (ER) and 
SF networks as the paradigm of homogeneous and heterogeneous networks respectively. By 
using FFS, we obtain lots of configurations at each interface along the nucleation pathway. 
From these configurations we scrutinize and compare the nucleating clusters in ER and SF 
networks. It is found that the processes of forming the critical nucleus are qualitatively 
different between the two cases of networks. For the former, a dominant cluster arise firstly, 
and groups smaller clusters gradually, while for the latter, many small clusters emerge at 
first and then abruptly turn into the critical nucleus. Interestingly, both the cluster size 
distributions follow power-law distributions and the slopes are nearly the same at early 
nucleation stage. 

The paper is organized as follows. Section [TTl presents the details of our simulation model 
and the numerical methods we employ to sampling the nucleation pathway. The numerical 
results are compared for SF networks and ER ones in Sec. IIIII A brief summary is given in 
Sec. Ml 

II. MODEL AND METHOD 

A. Network-organized Ising model 

We consider the Ising model on complex networks consisting of nodes. Each node is 
endowed with a spin variable Sj that can be either +1 (up) or —1 (down). The Hamiltonian 
of the system is given by 

if = - J ^ AijSiSj -h'^Si, (1) 

i<j i 

where J is the coupling constant and h is the external magnetic field. For convenience, we 
set J = 1 in the following discussions. The elements of the adjacency matrix of the network 
take Aij = 1 if nodes i and j are connected and Aij = otherwise. The degree, that is the 
number of neighboring nodes, of node i is defined as ki = J2f=i ^ij- 

The system evolves in time according to single-spin-fiip dynamics with Metropolis ac- 
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ceptance probabilities [30|, in which we attempt to flip each spin once, on average, during 
each Monte Carlo (MC) cycle. In each attempt, a randomly chosen spin is flipped with the 
probability min(l, e"''^'^), where /3 = l/^ksT) with ks being the Boltzmann constant and 
T the temperature, and Ai? is the energy change due to the flipping process. In the absence 
of an external magnetic fleld, the system undergoes an order-disorder phase transition at 
the critical temperature. Above the critical temperature, the system is disordered where 
up- and down-pointing spins are roughly equally abundant. Below the critical temperature, 
the system prefers to be in either of the two states: one state with predominantly up spins, 
and the other with almost down spins. In the presence of an external fleld, one of these 
two states becomes metastable, and if initiated predominantly in this metastable state, the 
system will remain for a signiflcantly long time before it undergoes a nucleation transition to 
the thermodynamically stable state. We are interested in the pathways for this transition. 

B. FFS method 

The FFS method has been successfully used to calculate rate constants and transition 
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29|, 
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32j. For 



paths for rare events in equilibrium and nonequilibrium systems 
clarity, we describe the method again here, together with some relevant details with our 
work. This method uses a series of interfaces in phase space between the initial and flnal 
states to force the system from the initial state A to the flnal state i? in a ratchetlike 
manner. Before the simulation begins, an reaction coordinate A is flrst deflned, such that 
the system is in state A if A < Aq and it is in state B ii \> Am. A series of nonintersecting 
interfaces A^ (0 < i < M) lie between states A and B, such that any path from Aio B must 
cross each interface without reaching Aj+i before Aj. The algorithm flrst runs a long-time 
simulation which gives an estimate of the flux escaping from the basin of A and generates 
a collection of conflgurations corresponding to crossings of interface Aq. The next step is to 
choose a conflguration from this collection at random and use it to initiate a trial run which 
is continued until it either reaches Ai or returns to Aq. If Ai is reached, the conflguration of 
the end point of the trial run is stored. This process is repeated, step by step, until Am is 



reached. For more detailed descriptions of the FFS method please see Ref. 33|. 

In this work, we will use FFS to study nucleation pathways of the equilibrium phase from 
the metastable spin phase. Speciflcally, we let > and start from an initial state with Si = 
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— 1 for most of the spins. We define the order parameter A as the total number of up spins 
in the network. The spacing between adjacent interfaces is fixed at 3 up spins. We perform 
1000 trials per interface for each FFS sampling, from which at least 200 configurations 
are saved in order to investigate the statistical properties along the nucleation pathway. 
The results are obtained by averaging over 10 independent FFS samplings and 50 different 
network realizations. 



III. NUMERICAL RESULTS 

^hose degree distribution 



34| . and the well-known 



In what follows, we employ a Barabasi-Albert SF network, w 
follows a power law P{k) ~ k^'^ with the scaling exponent 7 = 8 
ER random network. 

In Fig. [H we present schematically the evolution of local nucleating clusters in ER and SF 
networks at different stages (for clarity only show 100 nodes). Here, a nucleating cluster is 
defined as the component of connected nodes with up spins. Qualitatively, it shows distinct 
features along nucleation stages. In the ER case, there always exists a dominant cluster, 
which groups smaller ones gradually. While for SF networks, no dominant cluster appears 
at the early stage, but then a giant cluster emerges suddenly. This demonstrates that 
nucleation follows different pathways on ER and SF networks, indicating that heterogeneity 
of the network topology may play an important role. 

To further elucidate detailed characteristics along the nucleation pathway, we use FFS to 
generate configurations and perform detailed analysis on the nucleating clusters, including 
the largest cluster size, average degree of the cluster nodes, the number of clusters and cluster 
size distribution. According to CNT, there exists a critical nucleus size Ac of the new phase, 
above which the system grows rapidly to the new phase. Herein, we mainly focus on the 
nucleation stage where A < Ac. In our simulation, we determine Ac by computation of the 
committor probability Pb, which is the probability of reaching the thermodynamic stable 
state before returning to the metastable state. As commonly reported in the literature 
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2l|, the critical nucleus appears at Pb{K) = 0.5. Since Ac are different for different 
networks, we thus introduce A/Ac as the control parameter. 

For consistent comparison, we introduce S^ax as the ratio of the size of the largest 
nucleating cluster to over the total number of up spins, and plot Smax (averaged over the 
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FIG. 1: (Color online) Typical nucleating clusters for several values of A for the two different 
topologies studied (ER and SF). These shown networks contain 100 nodes, in order to have a 
sizeable picture of the system. Up spins and down spins are indicated by green circles and black 
circles, respectively. 

ensemble at each interface) as functions of A/Ac in Fig.2. Clearly, one can see that Smax 
for ER networks is always larger than that for SF ones. Specifically, at A/Ac = 0.5, S^ax is 
already more than 70% for ER networks, while it is only about 30% for SF ones, as shown 
by the dashed gray lines in Fig.2. But when A/Ac = 1 they almost tend to 100% together. 

To show our results more explicitly, we investigate the average degree Kn of the nodes 
inside the nucleating clusters, and plot Kn as functions of A/Ac in Fig. [3l As shown, Kn 
increases monotonically with A/ Ac for both ER and SF networks, which means the new 
phase tends to grow from those nodes with smaller degrees. For ER networks, Kn grows 
fast at the very beginning following by a relatively slow increasing. For SF networks, Kn 
increases slowly at first and jumps sharply when approaching the critical nucleus. Such a 
scenario is consistent with FiglHand FiglS 

To better understand the aforementioned differences, we present the number A^^ of the 
nucleating clusters as functions of A/Ac in Fig. Il](a). We observe that Ns non-monotonically 
depends on A/ Ac and the numbers of clusters in SF networks are always more than that 
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FIG. 2: (Color online) The relative size S^ax of the largest cluster as functions of A/A^. Parameters 
are N = 1000, the average network degree K = Q, h = 0.5, T/Tc = 0.3, Aq = 130 and Am = 880. 
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FIG. 3: (Color online) Average degree Kn of new phase nodes as functions of A/Ac. Other param- 
eters are the same as in Fig. 2. 

in ER ones. On the other hand, Ns for both networks approach the same magnitude near 
the formation of critical nucleus, but it decreases much more sharply in SF networks which 
is also consistent with the picture shown in Fig{T] to [31 In Fig. 111(b), the cluster size 
distributions P{s) for A/Ac = 0.2 and 0.3 are shown. Interestingly, P{s) follow apparent 
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FIG. 4: (Color online) (a) The number Ng of nucleating clusters as functions of A/Ac- (b) Size 
distribution P{s) of nucleating clusters, on a log-log scale, s denote the size of nucleating clusters. 
Other parameters are the same as in Fig. 2. 

power-law distributions in the small size range for both types of networks, and in addition, 
the exponents are nearly the same for fixed A/Ac- The power law breaks in the large size 
range, where large clusters dominate. 

The above results can be qualitatively understood in terms of CNT. CNT assumes that 
the formation of a nucleus lies in two competing factors: the bulk energy gain of creating 
a new up spin which favors the growth of the nucleus, and the interfacial energy cost, an 
opposing factor, which is due to the creation of new boundary links between up and down 
spins. That is, the change of free energy AF may be written as AF{X) = —2hX + aX, where 
a denotes the effective interfacial free energy, which mainly depends on the average number 
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of boundary links that an up-spin node has. Obviously, a node with more boundary links 
is more difficult to change its spin state. For SF networks, it is thus always easier for the 
leaf nodes with small degrees to change state than the hubs with large degrees. Since the 
degree distribution follows power-law, there exist a lot of hubs with intermediate degrees, 
as well as a few hubs with very large degrees. Usually, many leaf nodes are connected 
to relatively small hubs, which are further connected to large hubs. Therefore, only small 
nucleating clusters, consisted of leaf nodes and small hubs, can form at the early stage of the 
nucleation process. These small clusters are either away from each other on the network or 
separated by those crucial hubs with very large degrees. In the final stage of the nucleation, 
once the crucial hubs connecting these small clusters change their states, a giant nucleation 
cluster will emerge abruptly. This picture is consistent with those results shown in the above 
figures. For ER networks, however, the degree distribution follows Poisson distribution and 
no crucial hub exists, such that those new-formed clusters are usually connected together 
and one would not expect a sharp increase in the cluster size, which is observed in SF 
networks. 



IV. CONCLUSIONS 



In summary, we have studied nucleation pathways of the Ising model with Metropolis 
spin-fiip dynamics in ER and SF networks using the FFS method. Concerning the former, 
there always exists a dominant cluster which groups small clusters gradually until the critical 
nucleus is formed; while concerning the latter, many isolated small clusters grow separately 
which suddenly merge together into the critical nucleus. We have performed detailed analy- 
sis involving the nucleating clusters along the nucleation pathway, including the cluster size 
as well as its distribution, the mean degree inside the cluster, and so on, to further demon- 
strate the above scenario. The distinct nucleation pathways between ER and SF networks 
further emphasize the very role of network topology. Our study may provide a valuable 
understanding of how first-order phase transitions take place on complex networks, which 
could be of great importance not only for physical systems, but also for social and biological 
networks. 
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